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Abstract 

We present the micro-photoluminescence 
(jUPL) and micro-reflectance contrast (^RC) 
spectroscopy studies on thin Aims of MoSe2 
with layer thicknesses ranging from a mono- 
layer (IL) up to 5 L. The thickness dependent 
evolution of the ground and excited state exci- 
tonic transitions taking place at various points 
of the Brillouin zone is determined. Temper¬ 
ature activated energy shifts and linewidth 
hroadenings of the excitonic resonances in IL, 
2 L and 3 L flakes are accounted for hy using 
standard formalisms previously developed for 
semiconductors. A peculiar shape of the op¬ 
tical response of the ground state (A) exciton 
in monolayer MoSe2 is tentatively attributed 
to the appearance of Fano-type resonance. 
Rather trivial and clearly decaying PL spectra 
of monolayer MoSe2 with temperature co nf irm 
that the ground state exciton in this material 
is optically bright in contrast to a dark exciton 
ground state in monolayer WSe2. 

*To whom correspondence should be addressed 

^Laboratoire National des Champs Magnetiques Intenses 
CNRS-UJF-UPS-INSA, 25 rue des Martyrs, 38042 Greno¬ 
ble, France 

institute of Experimental Physics, Faculty of Physics, 
University of Warsaw, Pasteura 5, 02-093 Warszawa, Poland 

'ilCurrent address: Physikalisches Institut, Westfalische 
Wilhelms-Universitat Munster, Wilhelm-Klemm-StraBe 10, 
48149 Munster, Germany 


Introduction 

Layered semieondueting transition metal diehaleo- 
genides (TMDCs) sueh as M 0 S 2 , MoSe 2 , WS 2 , 
WSe 2 and MoTe 2 have beeome a subjeet of inten¬ 
sive study reeently due to their interesting optieal 
properties rieh in spin and valley physies, whieh 
largely results from the indireet to direet band gap 
erossover upon deereasing thiekness from bulk to 
the monolayer (lL).-“-^ Thin layers of these ma¬ 
terials have shown new interesting physieal phe¬ 
nomena sueh as valley polarization->^ and valley 
eoherenee^ effeets or tightly bound trions.-^^^ It 
has also been domenstrated that they ean serve as 
effieient sourees of light, either in the form of light 
emitting diodes— or single photon emitters, 
whieh is potentially promising towards teehno- 
logieal advanees sueh as in optieal eomputation 
and valleytronies. As sueh, the optieal properties 
of TMDCs are dominated by many exeitonie res- 
onanees over a rather small speetral region (a few 
eV). Therefore, a full eomprehensive pieture of the 
eleetronie band strueture requires a detailed exper¬ 
imental study of exeitons, and their exeited states 
as a funetion of layer thiekness and temperature, 
whieh has reeeived very limited attention so far. In 
this letter, we present a eomprehensive analysis of 
the evolution of many exeitonie ground state and 
exeited state features, spanning the pRC speetra of 
good quality IL to 5L thiek flakes of MoSe 2 . This 
information finds importanee for an improved the- 
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oretical understanding of excitons and charged ex- 
citons where large discrepancies have been found 
between the experimental and limited theoretical 
results (such as exciton binding energies).— 
Interestingly, for the A-exciton transition in our 
monolayer flakes, we find an unusual reflectance 
spectral lineshape, which is tentatively attributed 
to the appearance of a Fano-type resonance.—-— 
This could arise since the A-exciton lies in the en¬ 
ergy region spanned by the excited and the contin¬ 
uum states of the trion, where a strong interaction 
could exist.'^^*^! Temperature dependent evolution 
of the /iRC spectrum of the monolayer flake is 
consistent with the Fano hypothesis. Recent the¬ 
oretical works suggest that in addition to strong 
spin-orbit splitting in the valence band (VB) of 
TMDCs, there is also a splitting in the conduc¬ 
tion band (CB) (up to ~40 meV).— The sign 
of this splitting is expected to decide weather the 
energetically favorable exciton is optically active 
(bright) or inactive (dark) which has a strong in¬ 
fluence on the PL properties of the material. Our 
temperature-dependent /iPL studies on the mono- 
layer MoSe 2 are consistent with the theoretical 
predictions that the ground state exciton in this 
material is optically bright, which is found to be 
opposite to the case of WSe 2 monolayers.— We 
also derive the temperature activated energy shifts 
and linewidth broadenings of the ground state ex¬ 
citons from IL to 3L thick flakes and explain them 
on the basis of standard formalisms used for con¬ 
ventional semiconductors. 


Samples and characterization 

The thicknesses of the mono- to few-layer MoSe 2 
flakes were established using optical contrast, 
atomic force microscopy (AFM) and Raman spec¬ 
troscopy techniques. FigHJa) shows the Raman 
spectra obtained for the flakes under considera¬ 
tion. The Raman lines corresponding to the Ai^, 
E 2 g and B^^ phonon modes were observed in con¬ 
sistence with Ref. 0. The monolayer thickness 
was found to be ~0.65 nm using AFM (data not 
shown). The flakes were then characterized us¬ 
ing /iPL spectroscopy (see methods). Figure [U 
(b-d) displays the /iPL spectra of IL to 3L thick 
flakes. For the IL flake, two well separated fea- 



Figure 1: (a) Raman spectra obtained for the IL, 
2L, 3L, 5L and 180 nm (bulk) MoSe 2 flakes, (b) - 
(d) The low temperature /iPL spectra of the IL, 2L 
and 3L thick MoSe 2 flakes with the corresponding 
optical images displayed as insets. /iPL spectra for 
the 2L flakes appeared to be of two different types, 
shown by solid (optical image of the flake shown) 
and dashed lines (image not shown) in (c). 


tures corresponding to the emission from the n = 1 
ground state A-exciton and the charged exciton 
(trion, T) are observed. The emission intensity in 
the case of 2L flake is lower than for the IL by 
about 4000 times. A feature corresponding to the 
indirect gap emission (I) is observed, which is red 
shifted relative to the direct-gap emission peak A 
by 180 meV. Interestingly, in the /iPL spectra of 
some of the 2L flakes (solid line in Fig. [Itb)), we 
also observed two features which coincide with the 
A and T emission features of the monolayer, but 
are around 4 orders of magnitude weaker. These 
features existed throughout the area of rather large 
flakes (~30 jum wide), therefore the contribution 
due to the neighbouring IL flakes can be ruled out. 
Our hypothesis is that the process of mechanical 
exfoliation led to the formation of some pockets 
within the 2L flakes, which act as two separated 
monolayers giving rise to a monolayer like emis¬ 
sion. The juPL spectrum of the 3L flake was sim¬ 
ilar to that of the 2L flake where A and I emis¬ 
sion features were observed. The I feature was red 
shifted with respect to A by ~330 meV. We did not 
observe any measurable emission from the flakes 
with thicknesses larger than 3Ls. 
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Excitonic resonances versus 
layer thickness 

Figure [2l^a) shows the reflectance contrast spectra 
of the MoSe 2 flakes with layer thicknesses up to 
5Ls, except a 4L flake. The spectra are domi¬ 
nated by many excitonic resonance features. To 
determine the resonance energies, we fitted the 
spectra using a transfer matrix method which ac¬ 
counted for the interference effects on the exci- 
ton lineshapes due to the multilayered structure of 
the sample as reported earlier.— In our method, 
the excitonic contribution to the dielectric func¬ 
tion was considered to be a modified Lorentzian 
function— (see Methods section), where a phase 
factor /3 was used to account for the peculiar line- 
shape of the A-exciton feature in the case of IL 
flake (discussed later), and was set to be equal to 
0 for all other features (except A* for IL, where /3 
was assumed to be equal to that for A-exciton). It 
is tempting to derive the reflectivity and absorption 
spectra of the bare flakes (without substrate) from 
the spectral modeling. These spectra are shown in 
Fig. [Jtb) and (c) respectively where the absorp¬ 
tion coefficient a = 4nkjX is calculated by de¬ 
riving the imaginary part of the refractive index, 
k, from the modeling. One can clearly observe 
the features corresponding to the transitions tak¬ 
ing place at various points of the Brillouin zone, 
which have been identified earlier for bulk MoSe 2 
based on low temperature transmission measure¬ 
ments and theoretical calculations.— Following 
these reports, the observed transitions are named 
as A, B, A', B', d, C and D in the order of in¬ 
creasing energy. The shoulders A* and B* to¬ 
wards the higher energy sides of A and B exci- 
tons arise from the contributions due to the excited 
state exciton (n>2) transitions. An evolution of 
the transition energies with the layer thickness is 
shown in Fig.|2tb). The transition energies for the 
bulk MoSe 2 have been obtained from a previous 
report.— We notice that contrary to the case of 
WSe 2 ,“ where the A and B exciton resonances 
display a monotonic red shift when layer thick¬ 
ness is increased, their behavior in MoSe 2 is dif¬ 
ferent. The A-exciton feature shows an overall red 
shift (~58 meV) when the thickness is increased 
from a monolayer towards bulk, except in going 
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Figure 2: (a) Low temperature jURC spectra of 
the MoSe 2 flakes (solid lines) and the correspond¬ 
ing modeled curves (dashed grey lines) with thick¬ 
ness ranging from IL to 5Ls, except a 4L flake, 
(b) Reflectivity, and (c) the absorption spectra 
of the flakes derived from the spectral modeling 
in (a). The features representing A and B ex- 
citons, their excited states and the trion (T) are 
marked for IL flake, and other excitonic reso¬ 
nances are marked for 2L flake where they are 
more clearly visible. The spectra for thickness 
>2L have been shifted/amplified along y-axis for 
clarity by amounts shown above the curves on the 
left side, (d) shows the energies of corresponding 
exciton features as a function of layer thickness. 
Black diamonds represent the positions of broad 
features in the spectra, which possibly arise from 
overlaps between a few closely lying transitions. 


from 2L to 3L where it shows a slight blue shift 
(~5 meV). Furthermore, the B-exciton displays an 
initial blue shift from IL to 3L (~8 meV), which 
is then followed by a red shift (~13 meV) from 
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3L to bulk. Similarly, the A' and B' features show 
an initial blue shift from IL to 2L thiekness, and 
a red shift for thieker flakes. Although, a detailed 
theoretieal analysis is required to eompletely un¬ 
derstand this phenomenon, it ean be qualitatively 
explained by eonsidering the following two ef- 
feets: (i) a red shift of the band-to-band absorption 
edge and (ii) a deereasing exeitonie binding energy 
as the material undergoes a 2D to 3D transition, 
when the thiekness is inereased from a monolayer 
to the bulk limit.— The exeitonie transition ap¬ 
pears at an energy lower than the absorption edge 
by an amount equal to the exeitonie binding en¬ 
ergy. Therefore, depending upon the magnitudes 
of the two effeets, one ean either observe a red 
or a blue shift of the exeitonie transition. In our 
ease of the A-exeiton transition, the first effeet ap¬ 
pears to be mueh stronger than the seeond one re¬ 
sulting in a monotonous red shift. However, for 
B-exeiton (A' and B' exeitons), the seeond effeet 
dominates initially up to 3L (2L) thiekness result¬ 
ing in a blue shift, whereas the first one takes over 
later on whieh leads to a red shift. For the IL 
flake, a broad feature whieh appears around 2.7 eV 
(shown by diamond in Fig.l^b)) is most likely due 
to unresolved eontributions from d- and C- type 
exeiton features. A pair of d-exeiton features eould 
be resolved for thieker flakes whereas those related 
to the C-exeitons resulted in a broad feature around 
2.8 eV for 3L and 5L flakes (shown by diamonds 
in Fig. [2l^b)). D-exeiton features were identified 
around 3.2 eV for all the flakes but were not ae- 
eounted for in our modeling. Their positions were 
assumed to be at their minima and may be sub- 
jeeted to large errors. 

Now we diseuss the A* and B* features in Fig. [2] 
whieh we assoeiate with the speetral eontributions 
due to the n>2 exeited states of the A and B 
exeitons. The energies of these features agree 
quite well with those determined reeently using 2- 
photon absorption measurements.— Although the 
A* feature eould be observed for flakes with all 
thieknesses from IL to 5L, B* eould only be seen 
for the IL and 2L flakes and was too weak to be 
identified for thieker flakes. Also, only one fea¬ 
ture (2s-like) eould be identified for the A-exeiton 
for flakes of all thieknesses, whereas for the B- 
exeiton in the IL flake, both 2s and 3s features 
eould be observed, blue shifted by ~ 150 meV and 


~210 meV respeetively with respeet to the Is fea¬ 
ture. The exeiton series for IL B-exeiton signifi- 
eantly deviates from the 2D Rydberg model for a 
hydrogen atom in eonsistenee with reeent reports 
on IL WS 2 and WSe 2 ,— whieh renders it diffi- 
eult to prediet the exeiton’s binding energy with¬ 
out information about more (n>3) exeited state 
features. However, following the ideal 2D Ryd¬ 
berg model for the A-exeiton (eonsidering Is and 
2s states) and B-exeiton (eonsidering Is and 3s 
states), we find that their binding energies should 
be greater than ~100 meV and ~210 meV re¬ 
speetively. In analogy with our earlier findings for 
WSe 2 ,— the energy differenee between the 2s state 
feature shows an initial blue shift relative to the Is 
feature, in going from IL to 2L flake, for both A 
and B exeitons. This is followed by a red shift 
in the ease of the A-exeiton, when the thiekness 
is inereased from 2L to 5L. This phenomenon is 
explained on the basis of two eompeting effeets: 
(i) the non-loeal dieleetrie sereening of the exei¬ 
ton exeited states for the IL flake pushes the 2s 
state eloser to Is state,— and (ii) the material un¬ 
dergoes a 2D to 3D transition while transforming 
from the IL to the bulk form. The first effeet dom¬ 
inates when the thiekness is inereased from IL to 
2L whereas the seeond one takes over for thieker 
flakes resulting in a red shift. 

A-exciton in IL flake: possible 
evidence of Fano-resonance 

As illustrated in Fig. [Ifa to e), our standard anal¬ 
ysis of the refleetanee speetra (with /3 = 0°) pro¬ 
vides a eonsistent pieture of all observed exei¬ 
tonie transitions but elearly, this does not apply to 
the ease of the A-exeiton of the monolayer (with 
/3 ~ 90° at low temperatures, see below). The 
features related to all (but monolayer A-exeiton) 
imply the absorption resonanees in the form of 
Lorentzian peaks. In eontrast, the dedueed ab¬ 
sorption speetral shape of the IL A-exeiton res- 
onanee is unusual, displaying a dispersive-like 
form. It also shows a negative dip towards the 
higher energy side of resonanee [Fig. Ilte)], but 
for simplieity, we have not eonsidered any absorp¬ 
tion baekground in our modeling whieh will result 
in all positive absorption values. Among known 
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sources of unusual (non-Lorentzian) shapes of ex- 
eitonie resonances are the effeets of significant 
light-exciton coupling (polaritons) and/or the in- 
teraetion of the diserete resonanee with a eontin- 
uum of other states (Fano resonanee). Polaritonie 
effeets may be relevant in the refieetanee speetra 
of bulk materials— but ean generate only sub¬ 
tle effeets in two-dimensional systems—”^ unless 
the latter are embedded into a mieroeavity,-^^^ 
which is not the case of our struetures. We thus 
speculate that the observed unusual shape of the 
A-exeiton in IL MoSe 2 is due to Fano-type ef¬ 
fect. The exeitonie resonanees of the Fano shape 
ean indeed be observed in two-dimensional sys¬ 
tems. The best known examples are the Fano- 
type resonanees of the exeited exeitonie states eou- 
pled to the eontinuum of the ground state heavy- 
and light-hole exeitons in GaAs/GaAlAs quantum 
wells.—-— In search for a continuum of states 
which possibly couple to the A-exciton in the 
monolayer (and not in other layers), we note that 
it is only the monolayer whieh displays a promi¬ 
nent, below A-exciton absorption peak, related to 
the eharged exciton T. Obviously, there must ex¬ 
ist the exeited states of the T, ineluding the pre¬ 
dieted quasi-eontinuum of unbound T states whieh 
markedly eoincide in energy with the undressed 
neutral exeiton resonanee.— The effeetive eou- 
pling of the A-exeiton with the quasi-eontinuum 
of exeited states of the trion is our hypothesis to 
aeeount for the unusual shape of the A-exeiton 
in our IL strueture. Temperature dependent stud¬ 
ies eonfirm this eonjuneture on a phenomenolog- 
ieal ground. As shown in Fig. [3ta), the “unusual 
shape" of the IL A-exeiton transforms into a stan¬ 
dard form upon inereasing temperature in parallel 
to the progressive disappearanee of the eharged ex¬ 
eiton resonanee. This is aeeompanied with a grad¬ 
ual reduetion of /3 from ~90° to 0°. At tempera¬ 
tures above 200 K, when the trion is totally absent 
in the speetra (and the quasi continuum of its ex¬ 
cited states must have disappeared), the A exciton 
fully recovers its usual shape. 
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Figure 3: (a)-(e) Temperature evolution of the 
/iRC speetra in the speetral region eontaining 
A-exeiton features for IL, 2L and 3L MoSe 2 
flakes respeetively, along with the modeled eurves 
(dashed grey lines). The speetra for T > 20 K 
have been shifted downwards along y-axis for elar- 
ity. (d) The /3 parameter for A-exeiton in a IL 
as derived from the lineshape fits in (a), (e) The 
transition energy, and (f) the linewidths of the A- 
exeiton features as a function of temperature for 1, 
2 and 3 L flakes and A* feature in 3 L flake along 
with the fitting (solid lines) based on the models 
deseribed in the text. 

Temperature evolution of ab¬ 
sorption and emission spectra 

Fig. [3tb) and (c) show the temperature evolution 
of the A-exeiton /iRC features of the 2L and 3L 
thiek flakes. In the ease of a 2L flake at low tem¬ 
peratures (T<140 K), the main exeiton feature is 
aeeompanied with two small features at the higher 
energy side. Similarly to what has already been 
observed in the /iPL data (Fig.[TJb)), their energies 
eoineide with the A-exeiton and the trion in the IL 
flake. For the 3L flake, we also observe the evo¬ 
lution of the A* feature along with the A-exeiton 
in our speetral window. From modeling of these 
speetra, we deduce the transition energies and the 
linewidth broadenings of the exeitons as shown in 
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Fig-Hte) and (f) respectively. Similar to our results 
on WSe 2 ,— the temperature activated energy shifts 
can be equally well fitted using the Varshni— 
or O’Donnell et al— formulae. These relations 
are given by Eg{T) = Eq — {aT^)/{T + /3) and 
Eg{T) = Eq — S{hco)[coth{{hco)/2kT) — 1] respec¬ 
tively. In Varshni’s relation, a and /3 are the fitting 
parameters related to the temperature-dependent 
dilatation of the lattice and Debye’s tempera¬ 
ture respectively whereas in O’Donnell’s formula, 
{hco) and the S are the average phonon energy and 
coupling parameters respectively. These parame¬ 
ters for the A-exciton in the IL to 3L flakes and the 
A* feature in 3L are given in TablefUwhich behave 
very similarly to those we deduced for WSe 2 .— 
The linewidth broadenings of the A-exciton were 

Table 1: Fitting parameters as obtained from 
modeling of the resonance energies and lineshape 
broadenings for the IL, 2L and 3L ground state 
A-exciton features, and the A*(2s) feature in 3L 
flake. 


Parameter 

IL 

2L 

3L 

3L(A*) 

Varshni’s relation 

Eq (eV) 

1.658 

1.626 

1.632 

1.707 

a{l0-'^eV/K)5.61 

5.56 

4.80 

4.77 

/3(K) 

330 

330 

330 

330 

O’Donnell’s relation 

Eq (eV) 

1.656 

1.625 

1.631 

1.705 

{hco)imeV) 

19 

19 

19 

19 

S 

2.23 

2.18 

1.88 

1.88 

Rudin’s relation 

7o (meV) 

7.5 

12 

35.5 

45 

/ (meV) 

72 

93 

120 

92 


fitted using the relation proposed by Rudin et al.,— 
given as /(T) = 70 -f (JT 4 - //[l/(e^“/^^ — 1)]. 
Here, 70 is the broadening at 0 K, the term linear in 
T depicts the interaction of excitons with acoustic 
phonons (disregarded in present work being neg¬ 
ligibly small) and the last term arises from the in¬ 
teraction with LO (longitudinal optical) phonons. 
hco is the LO-phonon energy which was taken to 
be equal to 30 meV.- The fitting parameters for the 
four features are given in Table [B Similar to what 
has been observed for WSe 2 ,— / which quanti¬ 
fies the extent of LO-phonon assisted scattering of 
the ground state exciton with the excited states, in- 
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1.55 1.60 1.65 

Energy (eV) 


Figure 4: Temperature evolution of the /JPL spec¬ 
tra of the monolayer MoSe 2 . The spectra for 
T > 10 K have been shifted downwards along y- 
axis for clarity, and those for T>100 K are ampli¬ 
fied with the factors specified. The inset shows the 
integrated PL intensity for the exciton I^, trion ly, 
defect related emission Id, and the total emission 
intensity Iq as a function of temperature. 

creases from IL to 3L flake. This is because the 
number of excited states available for scattering 
increase with increasing flake thickness. However, 
/ is lower for A* than for A feature in the 3L flake. 
It suggests that the excited state feature has lesser 
number of available bound and continuum excited 
states for undergoing LO-phonon assisted scatter¬ 
ing, as compared to the ground state feature. 

Finally, we performed /iPL as a function of 
temperature from 5 K till 300 K. The results are 
shown in Fig. |4l The inset to the Fig. |4] shows 
the integrated emission intensity from the trion 
(If), the exciton (I^), and the total intensity (Iq, 
area over full spectral range) as a function of tem¬ 
perature. The remaining emission intensity [Iq- 
(I^-i-Ir)] is attributed to the defect related lumines¬ 
cence (Id), which contributes towards the low en¬ 
ergy side (below trion) of the spectrum and could 
result in a slightly asymmetrical trion emission 
lineshape.— The PL spectrum is initially domi- 
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nated by trion’s contribution at low temperature, 
which is quickly quenched with increasing tem¬ 
perature and becomes negligible when compared 
to the exciton’s contribution for T>200 K. For 
T>70 K, the exciton begins to dominate the spec¬ 
trum due to quick thermal redistribution of charge 
carriers in the k-space.— We notice that there is a 
slight increase in the lo (~Ia) from 200 K to 250 K 
(notice log-scale). Although this phenomenon is 
not well understood at this stage, we note that the 
change of the PL intensity upon increase of tem¬ 
perature results from the competition between the 
efficiencies of the radiative and non-radiative re¬ 
combination channels, the latter being particularly 
poorly known in the investigated structures. The 
observed rise of the PL intensity may therefore 
be a result of distinct temperature dependencies 
of these two recombination channels. Overall, Iq 
shows a net drop by a factor of ~ 400 when tem¬ 
perature was increased from 5 K to 300 K. This 
is in agreement with another recent report— and 
is opposite to what we have reported in the case 
of monolayer WSe 2 where Iq showed a gradual in¬ 
crease as a function of temperature.— It can be ex¬ 
plained as follows. It has been predicted theoreti¬ 
cally that the conduction band (CB) in the mono- 
layers of TMDCs also splits into two bands due 
to spin-orbit interaction.However, the sign of 
this splitting has been predicted to be opposite for 
MoSe 2 and WSe 2 . Therefore, the states of lower 
energy among these spin-orbit split CB states have 
opposite total angular momenta, due to which the 
lowest energy exciton transitions for MoSe 2 and 
WSe 2 are expected to be bright and dark respec¬ 
tively Typically, the integrated PL intensity in 
semiconductors shows a reduction when temper¬ 
ature is increased, due to increasing contribution 
from the non-radiative recombination processes. If 
the lowest lying CB state results in the formation 
of dark excitons in combination with the state at 
the top of the spin-orbit split valence band (VB), 
the PL intensity may be quenched at low tem¬ 
peratures where this radiatively inefficient state is 
preferentially filled with charge carriers. Upon 
increasing the temperature, the carriers may get 
transferred to the higher energy CB state upon 
thermalization, which opens up a radiatively fa¬ 
vorable channel for their recombination and may 
result in an enhancement of emission intensity 


(which has been observed for WSe 2 ). Therefore, 
contrary to WSe 2 , a gradual increase of jUPL inten¬ 
sity with increasing temperature is in agreement 
with the theoretical predictions that the lower en¬ 
ergy exciton in MoSe 2 is bright. 

Conclusions 

In conclusion, we studied MoSe 2 flakes with thick¬ 
nesses equal to a monolayer, 2Ls, 3Ls and 5Ls 
using /jPL and /iRC spectroscopy and as a func¬ 
tion of temperature. We observed the evolution of 
various ground-state resonance features called A, 
B, A', B', d, C and D as well as the excited-state 
counterparts of some of them i.e. A* and B*, as 
a function of layer thickness. The A-exciton fea¬ 
ture in the /iRC spectrum of the IL flake showed 
Fano-type lineshape at low temperatures. This hy¬ 
pothesis was supported using temperature depen¬ 
dence /iRC measurements. In addition, the tran¬ 
sition energies and linewidth broadenings of se¬ 
lected excitonic resonances in a IL to 3L thick 
flakes were determined as a function of temper¬ 
ature and their behavior was explained using the 
models routinely used for conventional semicon¬ 
ductors. Finally, our temperature dependent /iPL 
measurements are consistent with the theoretical 
predictions that the lowest energy exciton in IL 
MoSe 2 is bright. 

Methods 

The monolayer and a few layers MoSe 2 flakes 
were obtained on Si/(100 nm)Si02 substrate by 
polidimethylosiloxan-based exfoliation— of bulk 
crystals purchased from HQ Graphene. The flakes 
of interest were first identified by visual inspec¬ 
tion under an optical microscope. A transfer ma¬ 
trix based calculation was performed to predict 
the color of the flake as a function of its thick¬ 
ness.— The observed color of the flake was then 
compared with the theoretically predicted one to 
roughly estimate the number of layers in the flakes. 
Then the AFM characterization and Raman spec¬ 
troscopy were performed on the flakes of inter¬ 
est to unambiguously determine their thicknesses. 
The Raman spectroscopy was performed at room 
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temperature using an argon ion laser (514.5 nm, 
0.2 mW focused using a 50x microscope objec¬ 
tive), with a 0.5 m focal length monochromator 
and a 1800 lines/mm grating. For performing 
/rPL measurements, the 632.8 nm radiation from 
a He-Ne laser was focused on the flake using a 
50x long working distance objective. The sam¬ 
ple was mounted on the cold finger of a contin¬ 
uous flow liquid helium cryostat, at a temperature 
of ~5 K. The spot diameter was ~2 /im and the 
light power focused on the sample was 1 jUW for 
the monolayer, 500 /iW for 2L and 2 mW for the 
3L flake. The PL emission from the sample was 
dispersed using a 0.5 m focal length monochro¬ 
mator and detected using a liquid nitrogen cooled 
Si charge coupled device camera. For performing 
the /iRC measurements, the light from a 100 W 
tungsten halogen lamp was focused on a pinhole 
of 75 /im diameter. The light was then colli¬ 
mated and focused (spot size ~3 /im) on the sam¬ 
ple. The reflected light was detected using a setup 
similar to the one used for performing /iPL. If 
^(A) and .^o(^) are the wavelength dependent 
reflectance spectra of the MoSe 2 flake and of the 
Si/Si02 substrate respectively, then the percent¬ 
age reflectance contrast spectrum (A) was deter¬ 
mined using ^(A) = — ,^o(A)]/[^(A) -F 

,^o(A)]. The measurements were performed at 
temperatures ranging from 5 K to 300 K. For the 
lineshape analysis of the RC spectra, we followed 
a method similar to that described in Ref. [291 In 
this method, we considered the excitonic contribu¬ 
tion to the dielectric response function to be given 
by a modified Lorentz oscillator like model as 

£(E) = (nh + ik/,) + ^ ^2 _ £2 _ 


where n^, -|- iki, represents the background complex 
refractive index of MoSe 2 in the absence of exci- 
tons and was assumed to be equal to that of the 
bulk material— in the simulations. The index p 
stands for the type of exciton characterized by a 
resonance energy Ep, an amplitude Ap, a phase 
factor /3 (to account for the Fano-type lineshape) 
and a phenomenological broadening parameter jp 
[equal to full width at half maximum (FWHM) of 
the Lorentzian function]. The RC was then calcu¬ 
lated using the transfer matrix formalism. It may 


be mentioned that despite storing the samples in a 
vacuum desiccator, we noticed an overall drop in 
the PL intensity nearly by a factor of 3 when the 
PL measurements were repeated on the monolayer 
flake after 40 days which indicates the possible ef¬ 
fects of sample-aging. 
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